Microscopic model for multiple flux transitions in mesoscopic superconducting loops 
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A microscopic model is constructed which is able to describe multiple magnetic flux transitions 
as observed in recent ultra-low temperature tunnel experiments on an aluminum superconducting 
ring with normal metal - insulator - superconductor junctions [Phys. Rev. B 70, 064514 (2004)]. 
The unusual multiple flux quantum transitions are explained by the formation of metastable states 
with large vorticity. Essential in our description is the modification of the pairing potential and the 
superconducting density of states by a sub-critical value of the persistent current which modulates 
the measured tunnel current. We also speculate on the importance of the injected non-equilibrium 
quasiparticles on the stability of these metastable states. 
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When a superconducting loop is exposed to a perpendicular magnetic field the energy El of the L-th state is given 
by: 

where L — 1/2'k f \79ds is the vorticity, 9 is the coordinate-dependent phase of the superconducting order parameter, 
and the integration is made along the contour of the loop). S is the loop's circumference, $ is the magnetic fiux 
through the area of the loop, and 0o — h/2e is the superconducting flux quantum. If the system can relax to its 
ground state, then sweeping the magnetic field causes a periodic variation of the kinetic properties with period A$ — (f>Q 
corresponding to transitions AL = ±1. The persistent current in the loop is proportional to the derivative of the energy 
I ~ dEi^/d^, and shows the characteristic sawtooth behavior with the same period. The 'switching' supercurrent 
density corresponding to transitions from the state with vorticity L to the nearest state L ± 1 is j switch 1/S. This 
A$ = if)Q behavior is commonly observed at temperatures close to the critical onc"'^. 

Recently, low temperature experiments (T <C Tc) were performed on superconducting loops^*^*^ leading to jumps 
with vorticity changes of AL > 1. A phenomenological explanation, based on a numerical solution of the time- 
dependent Ginzburg-Landau equations, was given in Rcf.'^. The basic idea of that theory was as follows: at sufficiently 
low temperature and when sweeping the magnetic field up the system can be 'frozen' in a metastable state with vorticity 
Li. When the circulating persistent current reaches the corresponding critical density ~ 1/C this metastable state 
relaxes to the low-energy level L2 2> Li. Thus, in loops with perimeter S ^ ^ one may observe vorticity changes 
Ai ~ jc/j switch ~ 5'/^ 3> 1. One should notice that the explanation based on the Ginzburg - Landau formalism 
is strictly valid only at temperatures close to the critical one, and the extrapolation to the low-temperature limit 
requires further justification. Therefore, of particular interest are ultra-low temperature experiments^ where Al loops 
with circumference from a few /im to a few hundred iim were studied, and a well-defined periodic structure with very 
large vorticity changes AL up to ~ 50 were observed. In the present paper we will construct a microscopic model 
based on the Usadel equations^'^ to analyze the experiments of Ref.**. 

Within the Usadel formalism, superconducting correlations between electrons forming a Cooper pair are described 
by two complex functions: the pairing angle and the phase In the most general case, both functions depend 
on the coordinate r and the energy E. The pairing angle and the phase are linked through the coupled system of 
equations: 

^V^e+[iE- cos 9] sm9 + A cos 9 = 0, (2a) 

V(v, sin^ 61) = 0, (2b) 

where the superfluid velocity is Vs = DiVx ~ (2e/fi.)A), D is the diffusion coefficient and A is the vector potential. 
The relation between the normal state density of states at the Fermi level N{Q) and the superconducting density of 
states is given by N{E) = N {Q)Re[cos{9{E))\. The pairing potential A(r) should be defined self-consistently from 
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the integral equation: 

A(r) = 7V(0)Ve// / ta.nh{E/2kBT)Im{sme)dE, (3) 
Jq 

where Vef / is the pairing interaction strength and lod is the Debye frequency. These equations are vahd for diffusive 
superconductors with a mean free path £ much less than the 'pure limit' coherence length = O-lShvp/^o (where 
vf is the Fermi velocity of the electrons on the Fermi surface and Aq is the gap at zero temperature, zero magnetic 
field and in the absence of transport current). 

For loops consisting of a sufficiently narrow wire \/cr < ^ one may neglect the variation of the superconducting 
properties in the transverse direction, where a is the cross section of the wire. As a first approximation one may also 
neglect the influence of the tunnel contact on the properties of the superconducting loop'*. Within these approximations 
we obtain from Eqs. (2) and (3) a single algebraic equation: 

E + ircose = iA^^ (4) 
smt' 

where F = h{v'^)/2D is the depairing energy averaged across a. One may neglect the self-induced magnetic field (for 
the largest loop with S =100 ^m it leads to an error less than 1%), and consider the vector potential A as entirely 
due to the applied field A = {Hy/2, —Hx/2, 0). For arbitrary vorticity L in a rectangular loop the depairing energy 
F is: 

where the magnetic field H is measured in units of Hc2 = *i'o/27r^(0)^, w is the side of the loop and d is the width of 
the wire the loop is made of and both arc measured in units of ^(0) = y/HAo/D. 




0/(t)„ 

FIG. 1: Dependence of the gap (a), pairing potential A (a), current density (b) and jump in vorticity AL (inset in Fig. 1(b)) 
on the applied magnetic field in a superconducting loop with parameters: Tc = 1.2 K, ^(0)=150 nm, w=5 fim , d=120 nm, 
T=106 mK. 

It was shown^ that for quasi-one-dimensional superconducting loops of large circumference (S 3> ^) metastable 
states with fixed vorticity L„ become unstable if the persistent current density is equal to the depairing current 
density jc- In the present calculation we assume that when the current density in the loop (being a function of 
magnetic field H) reaches its maximum (at a fixed value of vorticity L) the loop switches to a new quantum state. 
The new value of the vorticity Lm is then found from the condition that the pairing potential A (minimum depairing 
energy F) for a given magnetic field is maximal. In terms of a quasiclassical description this transition criterion 
corresponds to the supervelocity Vs being equal to its critical value. In Fig. (1) the calculated dependencies of 
the superconducting gap, pairing potential A^ and superconducting current density are plotted as function of the 
applied magnetic field for a superconducting loop with parameters similar to the ones used in Ref^. Of particular 
interest is the fact that for small magnetic fields the period of oscillations is much larger than (po and this period 
decreases in higher fields. Both observations are qualitatively consistent with the recent experiment a^i^i^ . However, 
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simulations systematically give larger vorticity jumps than the ones observed in the experiment although provide the 
same dependence on the loop size (compare Eq. (6) and Fig. (7) in Refr4). Reasonable explanation of the discrepancy 
is the presence of inevitable imperfections acting as weak links in real samples. Even at the lowest temperatures the 
system exhibits transitions from the 'frozen' metastable states at smaller values of the magnetic field. In Ref.^ it was 
shown that an inhomogcncity in the loop can indeed strongly decrease the values of the actually observed vorticity 
changes. Close to the critical temperature thermal fluctuations disable the formation of metastable states, which 
results in the 'conventional' 0o periodicity (AL = ±l)J«Aft. 

It is interesting to compare the results of the present microscopic approach with calculations based on the Ginzburg- 
Landau model^i^. For a rectangular loop with side w and thickness d (w ^ ^ > d) at T/Tc <C 1 the transition criterion 
for a state with vorticity L = corresponds to the condition Fc/Aq ~ 0.235 (at T 0). Then the size of the vorticity 
jump is: 



AL ^ Nint W -r^— ~ Nint 0.68— 
I V Ao < / V < 



(6) 



where Nint(x) returns the integer value of the real variable x. The solution of the Ginzburg-Landau equations^ for 
the same transition gives: 



AL 



GL 



Nint 



1 2w 

71< 



Nint 0.58 



(7) 



The slight numerical discrepancy between Eqs. (6) and (7) is due to the different functional dependence of Fc in the 
Ginzburg-Landau model and the present one. The quantitative difference between Eqs. (6) and (7) is insignificant 
compared to experimental inaccuracies^iiiii. 
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FIG. 2: Dependencies of the tunnel current on the applied voltage: (a) at various strength of the tunnel barrier, zero magnetic 
field and fixed temperature, (b) at various magnetic fields, fixed temperature and tunnel barrier strength. 



So far we have considered persistent currents in an isolated superconducting loop. Periodic modulation of these 
currents by an external magnetic field can be measured using magnetization^''! or calorimetricS. methods. In the 
ultra-low temperature experiments of Refi^ an Al loop was used as the superconducting (S) electrode of a N-I-S 
junction being overlaped through a tunnel barrier (I) with a normal-metal (N) contact. The tunnel current Itun was 
measured at fixed bias voltage (or the voltage across the barrier at a fixed tunnel current) as a function of the applied 
perpendicular magnetic field. We argue that the observed periodic variation of the tunnel current (voltage) with 
magnetic field^ is due to the modulation of the gap and the pairing potential by the sub-critical persistent currentii 
in the loop-shaped superconducting electrode of the N-I-S structure. To calculate the tunnel current we use the 
conventional 'semiconductor' modeli^ which gives for the tunnel current: 

Itun = a\Tf X 1^ ^[/(i?,T) - f{E + eV,T)]dE (8) 
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where f{E,T) is the Fermi distribution function, a is the constant of proportionaUty, and T is the tunneUing matrix 
element. It can be shown that ckjTp ~ l/eRx, where is the tunnel resistance of the junction in the normal state 
and e is the electron chargei^. In our case of strong currents, to obtain the density of states in the superconductor 
N{E) we solved Eqs. (2-5) self-consistently. 

The described approach provides a qualitative description of the main features observed in the experiment of Ref.^^ 
multi-flux periodicity of the tunnel current and the decrease of the period of the tunnel current oscillation with 
increasing magnetic field. However, quantitative comparison is far from being perfect: the measured current of the 
voltage biased N-I-S junction in the limit eV < A is systematically higher than the value calculated using Eq. (8) 
(see Fig. 2(a)). We believe that to obtain a better quantitative agreement one should consider higher-order processes 
describing tunnelling of pairs of electronsiii^ giving rise to the sub-gap current. To account for this process we use 
a combination of the results obtained in Ref»i^ for pure superconductors (quasi-particle transport in the ballistic 
regime) and the simple 'semiconductor' modcl^^. This leads to the semi-quantitative result for the sub-gap current: 



a\T[ 



A{A, E)[f{E, T) - f{E + eV, T)]dE 



(9) 



where the probability of Andreev reflection A(A, E) is taken from Table II of Ref>i^: 



(E2 - A2)(^£;2/(-B2- A2) + 2Z2 + 1)2 



E <A, 
, E>A, 



(10) 



with Z being the dimensionless barrier strengtbi^, and A is found from solutions of Eqs. (2-5). The tunnel current 
Itun and the sub-gap current Igub have different functional dependencies on A, and hence their variation in magnetic 
field is also different. While sweeping the magnetic field the magnitude of the oscillations of the tunnel current is much 
higher than the magnitude of oscillations of the sub-gap current. Strictly speaking, this over-simplified introduction of 
a sub-gap current in our model is not applicable for dirty limit electrodes studied in Ref. 4. However, it demonstrates 
that at energies A << eV an account for the sub-gap current contribution gives qualitatively better agreement with 
experiment: Fig. 2 and Fig. 3. The described approach is not able to give quantitative exact values for the sub-gap 
current as other mechanisms (e.g. leakage current) might contribute to the total current measured in N-I-S structures. 
It is also known that the magnitude of the sub-gap current strongly depends on interference effects within the locus 
of the N-I-S junction, and hence is geometry-dependent^'^. 
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FIG. 3: Magnetic field dependence of the total current through a N-I-S junction for different values of the bias voltage (a) and 
normalized magnitudes of the low field {B < 5 mT) current oscillations Al/Imax as function of the normalized bias eV/A{T) 
(b). Parameters of the loop are the same as in Fig. 1. 



In spite of the obvious simplifications when describing the sub-gap current, Eqs. (8-10) combined with calculations 
of the density of states (Eqs. (2-5)), is able to give good qualitative and reasonable quantitative agreement with 
experiment-. As follows from Fig. 2, and in full agreement with the experimental findings^, the absolute value of the 
magnitude of the current oscillations as function of the magnetic field I{B, V = const) (Fig. 3(a)) and their normalized 
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magnitude AI / Imax{B,V = const) (Fig. 3(b)) depend on the bias voltage V. The only 'tuning' parameter used in 
fitting the calculations with experiment is the coefficient Z (see Fig. 2(a)). 

At high temperatures the current across a tunnel junction at <C A is mainly determined by the tunnel component 
Eq. (8) due to temperature smearing of the Fermi distribution function. While at very low temperatures the 
contribution of the tunnel current is negligible in the same limit, and the finite measured current is practically equal 
to the sub-gap term of Eqs. (9,10). For all temperatures and > A the total current is determined by the tunnel 
component. In view of the relatively weak dependence of the sub-gap current (see Eqs. (9,10)) on magnetic field, this 
explains qualitatively the existence of a maximum in the dependence of AI / Imax versus voltage at low temperatures, 
and its absence at higher temperatures (Fig. 3(b)). 

In order to improve our model one should include the injection of non-equilibrium quasiparticles from the normal 
electrode. Particularly at energies eV ^ A this effect results in a cooling of the normal-metal contact and a heating 
of the superconducting contact^^'-'^^. The inevitable consequence is a spatially- inhomogeneous modification of the 
superconducting gap within the locus of the N-I-S junction. Manifestation of such an effect has been observed in 
S-I-S'-I-S structures at temperatures comparable to the critical temperature of the S' electrod o^?i^? . At T — > there 
are very few equilibrium (thermally activated) excitations, and one might expect that even a small amount of extra 
quasiparticles injected from the normal electrode may affect the pairing potential A and, hence, influence the stability 
of the 'frozen' metastable states. We speculate that the mentioned non-equilibrium effect might explain also the 
variation of the period ABj of the current oscillation on the bias voltage (Fig. 6 and Fig. 8(c) in Rcf.^). The larger 
the applied voltage, the higher the quasiparticle injection, and, hence, the stronger the deviation of the distribution 
function in the superconducting electrode from its equilibrium value. Therefore a metastable state with a given 
vorticity becomes unstable at smaller values of the depairing energy F, or, in other words, at smaller values of the 
applied magnetic field. 

In conclusion, using a microscopic approach we analyzed the ultra-low temperature behavior of a mesoscopic-size 
superconducting ring in the presence of a magnetic field. The model was used to interpret recent experiments^ on N-I- 
S junctions with a loop-shaped superconducting electrode. The central result of the present work is the demonstration 
that the tunnel current oscillates in a magnetic field with a period, which scales with the loop circumference. For 
large loops flux changes are much larger than a single flux quantum. We found agreement with experiment. One 
should include the finite sub-gap current in order to improve the quantitative agreement between the present model 
and the measured total current across the junction. Using simple assumptions, we were able to obtain a reasonable 
quantitative agreement with the experimental results of Ref4.. 
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